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Abstract 

Using the supersymmetric quantum mechanics we investigate the wave 
function-sensitive properties of the supersymmetric potentials which have 
received a lot of attention in the literature recently. We show that a su- 
perdeep potential and its phase-equivalent shallow-partner potential give 
very similar rms vàlues for the weakly bound systems such as the deuteron 
and 11 Be nuclei. Although the corresponding eigenstates differ in the node- 
number, our investigation on the n Be(p,d) 10 Be single nucleon halo trans- 
fer reaction at 35 MeV show that also other physical quantities such as the 
cross section angular distributions calculated using these wave functions 
reflect the nodal structure rather weakly. This lends support to two nearly 
equivalent treatments of the Pauli principle. 

PubUshed m Eur. Phys. J. A 9, 19-28 (2000) 

1 Introduction 



The long-standing dichotomy of choosing a shallow or deep effective local poten- 
tial to describe nucieus-nucleus elàstic scattering was greatly clarified by Baye 
jU, who demonstrated that in scattering, for all practical purposes, these two 
kinds of potentials in the literature @],f| are phase-equivalent supersymmetric 
partner potentials of each other. Supersymmetric quantum mechanics allows one 
to transform a hamiltonian to its partner such that they possess identical spectra 
or differ at most from each other by having the lowest eigenstate eliminated in 
its partner's spectrum Q. Repeated application of the supersymmetric trans- 
form could yield a hamiltonian which has a prescribed number of eigenstates 
less than the starting hamiltonian. By the same token, a supersymmetric trans- 
form can also add an eigenstate of the desired energy to the starting spectrum. 
The supersymmetric transform, apart from eliminating the lowest eigenstate, in- 
duces in the first instance a change in the phase shifts of the continuum states 
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of the starting spectrum. Baye JI| showed that by a judicious choice of repeated 
supersymmetric transforms, the partner hamiltonian can be made to be phase- 
equivalent. Thus, in the case of a — a scattering, where the starting hamiltonian 
has a deep effective local potential H known to possess bound states which are 
unphysical because of the Pauli principle, Baye eliminated these unphysical states 
by successive supersymmetric transform while preserving the phase equivalence. 
In the course of this process, the original potential is transformed to be singular 
and progressively shallower. Moreover, the supersymmetric transforms are per- 
formed for each angular momentum (£) separately. Therefore, at the final stage, 
a set of shallow, energy-independent and ^-dependent effective local potentials 
is obtained which is completely phase-equivalent to the original deep, energy- 
independent and ^-independent potential. The set of shallow potentials bears a 
remarkable resemblance to the ^-dependent at — at potentials documented in the 
literature 0]. 

Such ambiguous choice in the nature of the potential to describe nucleus- 
nucleus elàstic scattering is found not only in a — a but also in many other 
systems 0. This dual picture of the nucleus-nucleus interaction arises most 
probably because the many-body description of the system has been simplified in 
different ways to a two-body interaction between two structureless partides. For 
the analysis of elàstic scattering, such dràstic difference in the character of the 
potentials is immaterial since only phase shifts are required and these potentials 
are phase-shift equivalent or in the main so. However, when one has to choose 
one of these potentials to be used in nuclear structure studies in which wave 
functions are explicitly involved, the implication of a deep or shallow potential 
may be immense. The supersymmetric procedure produces two sets of wave 
functions for weakly bound systems, which coincide at large distances but differ 
at small distances by the additional node appearing inside the core by use of 
the deep potential. Therefore, it is important to have a quantitative criteria to 
discriminate which potential has the correct wave function-sensitive properties. 

These deep and shallow potentials were used independently by Baye et al. and 
Liu || to calculate bremsstrahlung emission in possible nucleus-nucleus col·lisions. 
The calculations seemed to indicate that while the bremsstrahlung cross sections 
from resonating-group method and the deep potential resemble each other, those 
of the shallow potential are distinctly different. Hence it was concluded that the 
deep potential is to be preferred over its shallow partner if wave function-sensitive 
properties are important. 

Recently Dijk et al. |7j, and Ridikas and his co-workers JF|, have separately 
shown that a superdeep potential and its supersymmetric partner give very simi- 
lar rms vàlues for the model deuteron, and one-neutron halo systems (considering 
n Be nucleus), respectively. In the work of Ridikas et al. it was also stated that 
other physical quantities, which are more sensitive to the behavior of the radial 
wavefunctions in the nuclear interior, such as transition probabilities obtained by 
the deep potential and its phase equivalent partner, reflect the nodal structure 
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rather weakly. 

In order to complete our discussion, we repeat these calculations with some ex- 
tent, considering a transfer reaction, n Be(p, d) 10 Be, involving the weakly bound 
deuteron and single nucleon-halo 11 Be nuclei which are well suited for study- 
ing the consequence of different wave functions from the deep and reconstructed 
phase equivalent shallow potentials on reaction observables. Current experimen- 
tal activity in the area of light-neutron rich and drip-line nuclei now dictates the 
ràpid development of calculable theoretical models for reactions and scattering 
of effective few-body systems. Hence, there is an increasing general interest in 
supersymmetric potentials in this context. Our results are contributions to the 
discussion in this subject -an investigation in a relatively unexplored area of the 
quantitative consequence of the supersymmetric potentials. 

We begin with a brief sketch of the general method in section 2, where we 
also introduce a two-parameter, deep, shape invariant two-body potential used 
throughout the present calculations. In sections 3 and 4, we discuss the ap- 
plication of the method to the deuteron and 11 Be nuclei respectively, giving the 
characteristic properties of the constructed phase-equivalent two-body potentials, 
and the connection between the exclusion of deep-lying Pauli forbidden bound 
states from some potential and supersymmetry is reviewed in the light of the 
calculation results. Section 5 discusses the u Be(p,d) 10 Be reaction calculations 
in terms of the supersymmetric partner potentials. Finally section 6 contains a 
summary and the conclusion. 

2 Supersymmetric quantum mechanics 

Supersymmetric quantum mechanics || and its connection to the factorization 
method || have been extensively investigated JTÍJ. Since the ground state wave 
function v|/( ri=0 ) for a bound system has no nodes, it can be written as 



*("=°)(r) = exp(-^ I W(r)dr) (1) 
Introducing the operators 

B = W(r) + -!=p , B + = W(r) - -%==p (2) 



the hamiltonian can be easily factorized 

H - £ (n=0) = B + B (3) 

where E^ n=0 ^ is the ground state energy. Since the ground state wave function 
satisfies the condition 

B #( n =°)\ = (4) 
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the supersymmetric partner hamiltonians (H m , m = 1,2, ...) 

H\ = B B , H 2 = BB^ 



(5) 



have the same energy spectra except the ground state of Hi, which has no corre- 
sponding state in the spectra of H. The corresponding supersymmetric partner 
potentials are given by 



V^r) = [W(r)f 



h dW 



2/i dr 



V 2 (r) = [W(r)f + 



% dW 



2/i dr 



(6) 



It was shown that a subset of the potentials for which the Schròdinger equa- 
tions are exactly solvable share an integrability condition called shape invariance 
OH . The partner potentials of Eq. (6) are called shape invariant if they satisfy 



the condition 

V 2 {r-a l ) = V 1 {r-a 2 ) + R{a 1 ) (7) 

where ai j2 are a set of parameters that specify space-independent properties of 
the potentials (such as strength, range, and diffuseness) , a 2 is a function of a\, 
and the remainder R(ai) is independent of r. 

An iterative procedure within the supersymmetric quantum mechanics frame- 
work for building the partner of a given potential admitting the same eigenvalues 
except for that of the missing ground state was proposed by Baye |] on the basis 
of a general procedure due to Sukumar ||, ]Ï2[. The method relies on a factor- 
ization property of the hamiltonian, and makes possible the (exact) construction 
of the partner potential starting from the original potential ground state wave 
function. It actually requires two steps, the intermediate potential (V2) having 
the same negative energy spectrum for a bound system as the original potential 
(Vi), except for the ground state of the latter, but a different phase shift; V 2 {r) 
is given by 

h 2 d 2 

V 2 (r) = VAr)-2 — —\nm 1 



2/irfr 2 



where ^^E^' 



jx ) denotes the original ground state wave function. 
step provides V%(r) the final phase-equivalent potential (PEP) as 



(8) 



The second 



V3M = Vy(r) 



2/i dr 2 



^ 1 (E[ n - 0) )^ 2 (E[ 



(n=0)^ 



(9) 



where ty 2 (E[ n °') stands for the wave function at the same energy, but calculated 
with the intermediate potential V 2 . Eq. (9) can also be reduced to the form 

h 2 d 2 



V 3 (r) = V^r) - 2- 



ln 



(n=0) 



dr' 



(10) 



2/i dr 2 

Elimination of more than one state is accomplished by iterating this two-step 
procedure. 
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2.1 A two-parameter superdeep potential 

In Ref. 0], the on-shell equivalence of the deep quantum-chromodynamically 
motivated realistic nucleon-nucleon interaction proposed by Kukulin et al. HÏ4| 
with more conventional repulsive-core forces has been investigated by eliminating 
its unphysical deeply bound states, while preserving its scattering properties and 
the binding energy of the deuteron. Using the spirit of this work, and of Ref.0, 
here we use an alternative superdeep potential. As a simple, physically interesting 
example, consider the potential 



V(r) = —Vosech 2 (5r 
Potentials of this shape can be generated from the superpotential JÏ5 

W(r) = Atanhpr, A > 0. 
In fact, using Eq. (6), the supersymmetric partner potentials are 

/3h ' 



(12) 



Vx(r;A) = A 2 -A\^A + 
V 2 (r; A) = A 2 - A \A - 
Clearly, one can write 



sech 2 (3r , 
| sech 2 (3r . 



(13) 



V 2 {r-A) = V 1 {r-A 



(3% 



+ A 2 - A [A 



2 fi 



(14) 



which is precisely the requirement of Eq. (7) for shape invariance. Therefore the 
bound state energies of the potential V\ are 



El 



(«) 



À 



A-n 



2n t 



(15) 



The energy levels £ (n) of the original potential given by Eq. (11) can be 
obtained by subtracting A 2 from and identifying 



V = A[A + 



Solving for A and requiring A > gives 



ph 



2u 



(16) 



.4 



ph 



_ pn?_ 

'2y/2jl ^ 2\ ^2JI 



1 

+ ^ 



+ 4VÍ, • 



(17) 
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Therefore, the energy levels of the deep potential V(r) = —V^sech 2 (3r are 



E {n) = E H _ A 2 = _( A _ 
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which is well known to be the correct answer ||16| . As we deal with the bound sys- 
tems, we require the odd solutions due to boundary conditions. Hence replacing 
n in Eq. (18) by 2n + 1 term, we arrive at , 

E {n) = -^-{À-2n- l) 2 3 2 , n = 0, 1, 2, ... , (19) 

where A = j0rj^ ■ The depth of the potential given by Eq. (16) reduces in this 
case to the form 

V = -^À(À + 1)3 2 (20) 
2// 

Eqs. (19) and (20) are in consistent with the expressions used in Ref. [[7] 
where also the analytical expressions for the wave functions of the ground and 
first excited state corresponding the potential of interest can be found. 

We have first employed this deep sech-squared potential in analyzing the 
alpha-alpha scattering (but not discussed here) by choosing the two parame- 
ters as À = 5.945 and 3 = 0.535 frrT 1 , together with = 10.375 MeV fm 2 . 

We have observed that -V Q sech 2 3r -> -Uoexp^^ with U = 122.694 MeV, 
a = 0.22 fm~ 2 and have reproduced successfully the Figs. (1,2) of Ref. |T[] using 
this two-parameter shape invariant superdeep potential, without involving a gaus- 
sian type potential in the calculations. However, as we deal with u Be(p,d) 10 Be 
reaction calculations throughout the present work, we focus on the treatment of 
deuteron and u Be ground state wave functions in the supersymmetric quantum 
mechanical framework. 



3 Application to the deuteron system 

All the available "realistic" nucleon-nucleon forces are characterized by a rela- 
tively weak central attractive part and by the presence of a hard or soft repulsive 
core at small distances; the first feature reflects the loose binding of the neutron- 
proton system, while the introduction of a repulsive core is required by the nega- 
tive vàlues assumed by the low-experimental phase-shifts when energy increases. 
However the feasibility of a description of comparable quality of the low energy 
properties of the two-nucleon system in the 1 Sq and 3 S\ - 3 Di channels (includ- 
ing deuteron properties), in terms of a deep, purely attractive interaction called 
Moscow potential was demonstrated by Kukulin and his co-workers |Ï4] |. Their 



potential differs from those obtained in the more traditional approaches by the 
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existence of an additional deeply-bound state in each channel, and by an increase 
of the absolute singlet and triplet phase-shifts due to this extra unphysical bound 
state. It is well known from cluster nuclear physics that these seemingly contra- 
dictory features - that is, repulsive core versus deep potential descriptions - are 
two ways to simulate the effects of the Pauli principle in a local potential model 
description when the two interacting partides are composed of identical fermions. 
There have been quite a few attempts to derive the features of the two-nucleon 
interaction from a quark picture of the nucleon. The most non-relativistic quark 
model calculations led to an effective nucleon-nucleon interaction with a strong 
repulsive core and an intermediate range attraction similar to those displayed by 
the empirical potentials. On the other hand, the work described in Ref. 0] in- 
dicated, for the nucleon-nucleon scattering, that the relative s-wave function has 
to have a node at small distance. The existence of a node in the relative motion 
wave function can readily incorporated in a local potential description, provided 
interaction is deep enough to accommodate one (nodeless) deeply bound state, 
such as the one proposed by Kukulin and his co-workers. The work described 
here will, in addition to the other investigat ions undertaken, demonstrate ex- 
plicitly the equivalence of such deep potentials with the more orthodox repulsive 
core empirical interactions, by constructing the phase-equivalent supersymmetric 
partner of the deep potential freed from unphysical bound states but still binding 
the deuteron with correct energy. The resulting supersymmetric potential for the 
deuteron case are to be shown to have a short range repulsive core followed by 
a shallow attractive part, which are very similar to those displayed by realistic 
interactions such as the Reid soft core potential Jïj 



3.1 Phase-equivalent potentials for the deuteron 

Now, superdeep potentials such as the Moscow potential give deuteron wave 
functions with a node at short distances. The node arises because there is an 
additional bound state which is Pauli forbidden for the actual neutron-proton 
{n — p) system. The latest version of Moscow potential |L4 1 includes both central 



and a tensor component, together with the central and tensor one-pion exchange 



contributions (OPEP). In addition, the equations in Ref.[14] to solve both for the 
bound state (S- and D-wave) and the scattering problem in case of triplet po- 
tential are more complicated, one needs to work out the usual coupled equations. 
For simplicity, in present calculations we use an alternative superdeep potential 
discussed in Section 2.1, which produces a model deuteron wave function that 
has a node like the Moscow potential. However, this simple potential does not 
have the required OPEP tail and do not consider D-wave related to the tensor 
potential. So the physical observables calculated by this potential, such as the 
radius of the deuteron which we will deal with later in this section, should not 
be compared with the experimental value. We note at this point that our aim 
here is not the rigorous reproduction of experimental data but to test the reli- 
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ability of a deep potential description involving some unphysical bound states, 
and in particular to search for the wave function sensitivity features of the deep 
and of its supersymmetric partner-shallow potentials used both in the entrance 
and exit channels of the u Be(p,d) 10 Be reaction. For this reason, the use of an 
appropriate simple potential, such as the binary sech-squared potential, in the 
present analysis does not cause any physical problem. 
Considering the well-known charge radius formula 

R\charge)= l -Rl + - A R 2 Tms (21) 

with R p being the proton radius and R rms the mean-square intercluster distance 
(matter radius), one can determine the free parameters A and f3 for the potential 
considered in analyzing the deuteron nucleus by solving the following system 
equations ; 

k 2 

#(») = (A — 2n — l) — -2.226 MeV , (22) 

2/^n— p 

Rrms = \jdrr 2 ^\À,{3,r) 2 = (1.95 fm) , 

where n denotes the energy level as stated earlier and m refers to m th hamiltonian. 
Throughout our calculations, ^ — is set 41.47 MeV fm 2 and the arbitrary 

constants A and (3, for the Moscow-type binary potential, are calculated as 3.146 
and 1.587 fm' 1 respectively. In this case, from Eq. (22), the ground state 
has a binding energy of about 481 MeV, which is unphysical and needs to be 
suppressed. The physically meaningful deuteron bound state for this superdeep 
potential corresponds to the first excited state having a binding energy of 2.226 
MeV. 

The building the partner of a given potential admitting the same eigenvalues 
except for that of the missing ground state has been discussed in Section 2. The 
two-parameter superdeep potential- V\{r) and its phase-equivalent supersymmet- 
ric partner-V^r), together with the intermediate non-phase-equivalent potential 
V^r), are shown in Fig. 1-a, and their corresponding wave functions in Fig. 2. As 
a result of the presence of one spurious bound state, the deuteron wave function 
for the superdeep potential possesses one node near the origin (around 0.56 fm). 
Figure lb compares the phase equivalent repulsive core interaction (V3) with the 
central part of the usual realistic Reid soft core potential JET]. The general fu- 
tures of the PEP (such as the radius of the repulsive core and the strength of 
the attractive part) are seen to be similar to those of the Reid soft core poten- 
tial. In spite of different analytical behavior near the origin for both Reid Soft 
core interaction (behaves at r — > as e -( constxx ) / x ) and the transformed phase 
equivalent shallow potential (V3(r) r ^ ~ (const/r 2 )), we observe in the figure 
a considerable similarity of both interactions. This means we have very tight 
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interrelation between a deep nucleon-nucleon model potential and the Standard 
Reid soft core interaction, which will be discussed later in this section. 

It is seen from fig. 2 that our reconstructed phase-equivalent supersymmetric 
partner potential (V3) has led to relative motion wave functions very similar to 
this generated by the deep potential ( Vi) outside the core region, but which lack 
the small distance radial node resulting from the suppression of the unphysical 
bound state. If there is a node in the wave function and the wave function is 
reasonably large at small distances, then one might expect that because of the 
normalization the wave function at large distances would be reduced. In other 
words, the asymptotic phase-equivalent wave function will have a smaller value of 
the asymptotic amplitude and hence the radius would be reduced. To clarify that 
if these wave functions having different behavior inside the core lead to quanti- 
tatively different results, we investigate the dependence of radius calculations, as 
an observable, on the wave function properties. 

The deuteron matter radius R rms can be calculated numerically from Eq. (22) 
for either *(r) = V { " =1 \r) or *(r) = ^ n=0) (r), the bound state wave functions 
for the superdeep and transformed partner potentials. The numerical calculations 
of rms value (root-mean-square radius) show that the radius of deuteron is 1.953 
fm and 1.955 fm for the superdeep potential and PEP respectively, which are 
so close. That means as a physical quantity, the radius calculated using the wave 
functions having one-/no-node reflect the nodal structure rather weakly. But, in 
case of experiment requiring a reduction in R rms for a nucleus then short range 
contribution to the potential appears, such as the superdeep potential used here, 
would be necessary. 

In order to check the accuracy of the supersymmetric quantum mechanical 
methods used in constructing the phase-equivalent potential, we have carried 
out additional calculations on the phase-shift. The resulting phase-shift curve 
obtained by PEP (Vj(r)) is compared with that obtained by the deep (Vi) po- 
tential and an excel·lent agreement between the resulting and initial phase-shifts 
is observed. 

To summarize, when nucleons are endowed with quark structures, nucleon- 
nucleon bound and scattering properties can be described by a deep potential 
(ps 1000 MeV), whose supersymmetric partner potential is singular and looks 
teasingly like the Reid soft-core or similar shallow potentials. One can make a 
conclusion at this stage on the supersymmetry aspect of nucleon-nucleon interac- 
tion in the light of the calculation results obtained, making a possible connection 
between the exclusion of deep-lying bound states from some potential and super- 
symmetry. This very interesting aspect is connected with very deep interrelation 
between many-body and potential treatment for the composite particle interac- 
tion in case of the partides composed of elementary fermions. From physical 
point of view this interrelation means the existence of a deeply hidden relation 
between the relative motion of composite nucleus comprised from fermions and 
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the internal excitations of the composites. In fact, in non-relativistic quantum 
mechanics which is used for the treatment of interaction of the composite par- 
tides, the relative motion of the composites is treated as a bosonic degree of 
freedom (i.e., no Pauli principle constrains are put to the relative motion). On 
the other hand, the internal excitations of the quark degrees of freedom inside 
the composites in the process of the mutual collision of the composites should be 
treated as the manifestation of fermionic degrees of freedom. The main problem 
in description of composite particle interaction is the complicated interrelation 
between relative motion of the composites and their inner excitations. And from 
this point of view, the existence of the above supersymmetry aspect could mean 
that the collision of such composites should be described correctly only within the 
framework of supersymmetrical quantum mechanics and is only the projection of 
this nontrivial picture onto mutual relative motion of the composites. 

As a conclusion, by means of supersymmetric quantum mechanics we have 
found that the two-parameter superdeep potential and its supersymmetric part- 
ner, which is phase-equivalent to the former and freed from the unphysical deeply 
bound states, give very similar rms vàlues for the deuteron system. Although 
the corresponding eigenstates differ in the node-number, our investigations have 
shown that the matter radius calculations using these wave functions reflect the 
nodal structure rather weakly. This lends support to two nearly equivalent treat- 
ment s of the Pauli principle by choosing the physical solution either by node- 
number criteria or by inclusion of a repulsive part of the potential at the origin. 
However, the similar rms results obtained do not automatically imply that other 
observables such us differential cross section, vector analyzing power calculated 
using these wave functions have to coincide. Clearly considerable additional work 
is still needed to test further the virtues of a deep potential description of the 
nucleon-nucleon interaction, which will be discussed in detail in Section 5 con- 
sidering a halo transfer reaction. However, as the ll Be nucleus in the entrance 
channel of the reaction considered is weakly bound, like the deuteron in the exit 
channel, a similar discussion for the 11 Be system being one-neutron halo nucleus 
within the framework of supersymmetric quantum mechanics is necessary before 
proceeding. 

4 Application to the n Be system 

Research with radioactive nuclear beams is currently one of the most active areas 
in nuclear physics. As one of the successful applications of such nuclear beams, 
exòtic structures have been observed in nuclei near to the driplines, which are 
called, as the most interesting discoveries, the neutron halo. These nuclei have 
opened studies of weakly bound nuclear systems, which has not freely accessi- 
ble before. These nuclei, such as ll Be, have long-range wave functions and are 
characterized by a cloud, or halo, of neutron probability that extends far beyond 
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the dense core. According to classical physics such nuclei should not exist at 
all because the strong nuclear force (the glue that binds neutrons and protons 
together) has too short to hold the far off neutrons in the halo. Instead, they 
owe their existence to quantum theory which describes the location of subatòmic 
partides by a mathematical cloud of probability. In [jl8| and references therein 
essential features of loosely bound systems, having an unusually large size, are 
discussed. In this Section, as an example to halo systems, we examine the ground 
state of u Be, which consists of a single neutron halo with a 10 Be core nucleus, 
using the supersymmetric quantum mechanics. It is well known that the domi- 
nant component of the n Be ground state is produced by the coupling of a lsi/2 
neutron, having the separation energy of 0.503 MeV, to a 10 Be core. 

4.1 Phase- equivalent potentials for the 11 Be system 

In this analysis, which will be restricted to the s-motion only, we again make use 
of the sech-squared deep potential with appropriate physical parameters. At this 
stage we should stress that one can relate the matter root-mean-square radius 
RRMs{ ma tt er ) to the single neutron root-mean-square radius R rm s by the formula 
^ 

w w 

R RM s( matter ) = Jyy + ^ R RMs( COre ) + f W+1 \2 R r™ ^ 

where W is the mass number of the core, here the l0 Be nucleus. Using the 
RnMs( core ) — 2.3 fm and RRMs( ma >tt er ) — 2.73 fm, it is easy to check that 
the value R rms = 6.70 fm gives roughly the average of the measured vàlues of 
Rrms (matter). 

To determine the two free parameters (A, j3 ), one needs to solve the following 
equations leading to the correct rms value and binding energy for the Be — 11 
system, 

h 2 

E (n) = (A — 2n — l) = —0.503 MeV; (24) 

2 fi 

Rrms = J drr 2 ¥^(À,P,r) 2 = (6.70 fm) 

where /i is the reduced mass of the system: 10 Be + n. 

The physical solution is chosen using the node-number and parity quantum 
number criteria. As we take into account that Os 1/2 the orbit is completely 
occupied and put the single neutron in lsi/2 state (having one node in the wave 
function) to have a positive parity required by the experiment, we choose the 
excited state wave function with one node, instead of ground state, as a physically 
meaningful solution for the superdeep potential in analyzing the ground state of 
11 Be nucleus. 
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Now we are ready to apply the supersymmetric technique to create the PEP 
corresponding to the sech-squared superdeep potential for the 11 Be one-nucleon 
halo system, with the calculated vàlues of A — 3.124, and (3 = 0.694 jmT 1 to- 
gether with ^ = 22.81 MeV fm 2 . Using the formulae given by the previous 
sections, the supersymmetric partner potentials and corresponding eigenstates 
are calculated. Fig. 3 illustrates the superdeep and related partner potentials 
while Fig. 4 gives their eigenstates respectively. It is seen from Fig. 4 that the 
eigenfunction corresponding to the PEP is a nodeless ground state. In this case 
the Pauli principle is taken into account by the repulsive part of the potential, 
see Fig. 3, repulsive up to 1.5 fm approximately. 

If one calculates the R rms for the system of interest using the nodeless eigen- 
state of PEP, the value of 6.78 fm is obtained. Even the wave functions differ 
in the node number, like the deuteron case discussed in Section 3, the R rm s is 
nearly the same for the phase-equivalent potentials, which is found about 6.70 
fm for the initial deep potential, whereas 6.17 fm for the intermediate non- 
phase-equivalent potential. One should not forget that the non-PEP potential 
does not belong to the PEP family, that is why the corresponding R rms vàlues so 
close for the PEP potentials while that of non-PEP considerably differ than the 
others. Moreover, the corresponding eigenfunction of this non-PEP potential has 
a different asymptotic behavior as well. 

We have also calculated the s-wave phase shifts in case of 10 Be(d,p) ll Be 
elàstic scattering up to 20 MeV for the three potentials. The results clearly 
exhibit the difference between PEP and non-PEP. 



5 Application to the n Be(p,d) w Be reaction at 
35 MeV 

Single nucleon transfer reactions, such as the (d, p) and (p, d) reactions, have been 
a reliable tool in nuclear spectroscopic studies of stable nuclei, determining posi- 
tions, spins and parities of nuclear states. Recently, the use of low energy single 
nucleon transfer reactions for structure studies of exòtic nuclei have attracted 
attention [r9]-||20[|. Because of the simplicity of the theoretical interpretation of 



these reactions, they are thought to provide an important source of the infor- 
mation about the structure of halo nuclei, such as 11 Be. It is now understood 
that the lsi/2 neutron single particle state in this region is lowered and that a 
dominant component of the ll Be ground state is produced by the coupling of a 
lsj/2 neutron to a 10 Be (g.s., + ) core; with a smaller but significant component 
in which a neutron is coupled to a 2 + excitation of the 10 Be core. 

The importance upon such transfer reaction spectroscopic studies of the in- 
clusion of the deuteron breakup degrees of freedom has been well discussed in 
Ref. EO] via the theories used to analyze measured cross section observables, and 
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shown that the magnitudes of the calculated cross sections, and particularly the 
ratio of the cross sections to the ground state and 2 + core states, of 10 Be are 
affected by the inclusion of three-body channels. 

Here in this section we do not discuss the details of these calculations. The 
aim of the present calculations is to investigate how the calculated physical ob- 
servables of the reaction sensitive to the nodal structure, involving the deuteron 
and bound neutron wave functions discussed in the previous sections with one 
node/no node. 

We calculate the transfer amplitude using the prior form of the (p, d) matrix 
element, thus the transition interaction is the n — p interaction and we need a 
full (three-body) description of the n+p+ 10 Be system in the final state. For the 
description of this final state we have used both the adiabatic (AD) model pl| 



and the quasi-adiabatic (QAD) approach However, for the sake of clarity 
in discussing the physics behind this application, we here consider only the AD 
model calculations as both model calculation results have led us to the same 
conclusion. 

In the context of (p, d) reactions the outgoing deuteron wave function, \E^ ; 
enters the transition amplitude 

T p d = (^d(r, R) V np \xp(r p )(f)n(rn)} , (25) 

where r (= f p — r n ) is the relative coordinate of the n—p pair and R — | (f p + r n ) 
is the center of mass coordinate. Here <p n is the neutron bound state and Xp 
incoming proton wave function. The vectors r n and r p are the positions of the 
transferred neutron with respect to the l0 Be core and of the proton relative to 
the ll Be system. The cross section (in the center of mass frame) for neutron 
pickup to bound n — p pair is then given by 



do v d _ HiHf kf 



T pd \ 2 (26) 



where dílf is the element of the sòlid angle for the asymptotic center of mass 
momentum of the bound deuteron. In the above pi and pf are the center of mass 
relative motion reduced masses in the initial and final channels, ki and kf are the 
entrance and exit channel wave numbers, respectively. Throughout this paper we 
restrict the formalism to s-wave n—p relative motion for simplicity. In zero-range 
approximation then it is the wave function at coincidence, ^ (r « 0, Rj which 
is of import ance. 

5.1 Transition amplitudes in zero-range approximation 

It is important for the present work to make clear the essential differences in 
the calculations carried out using both phase-equivalent potentials, the initial 
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superdeep sec/i-squared potential and the shallow phase-equivalent partner po- 
tential, in calculating the ground state deuteron wave function in the final state 
and for the calculation of the transferred neutron bound state wave function in 
the initial channel of the reaction. So they are developed in some detail within 
a common notation. The calculations employed the deep potential (denoted by 
V\ in sections through 2 — 4) will be represented by the script deep in the follow- 
ing formulae and for the phase-equivalent shallow potential (represented by V 3 
earlier) calculations we use the script pep. 

The transition amplitudes for the processes are evaluated in a zero-range 
approximation, and the related relative amplitudes, of primary interest here, 
are accurately described. For clarity we will not show the transferred neutron 
spectroscopic factor or any spin projection labels explicitly. 

In the adiabatic approximation the required transition amplitudes are 

d ee PÍpep) = / Xd ( r ,R)$f e ^\r)\ YS^\r) \x P (r P )<i> d n eep{pep \r n )) (27) 



As in calculating the transfer amplitudes, we make use of the zero-range 
approximation, thus for the bound deuteron we replace 

V^^\r) Xd (r,R)4 ee P(peP \r) « D d eepipep \p, d) Xd (r » 0,R)5(r) (28) 
with strength parameters 

D d eep(pep \p, d) = V^J drrV^ ep \r)ut ep{pep \r) , (29) 

where uo is the radial deuteron ground state wave function (£ np = 0) such that 
$d(r) = The similar replacement should also be done for the QAD 

calculations. 

It follows that, 

T dee P (pe P ) = pdeepipep) ^ Q M deep(pep) fa q ^ 

= Dt P{pep \p,d)( Xd (r*0,R) Xp (^ n ^\R)) 



where 7 = t^tt- Finally, considering the general expression for the cross section, 



we obtain 



7 dee P ( P e P ) , 

d(T P d _ IW± k f 2 



D 2 (deep,pep) \M deep{pep) (p, d) (31) 



dn f (2-nh 2 ) h 

for the calculations used the adiabatic model in describing the final state. 
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5.2 Calculation methods 



We calculate the cross section angular distributions for the n Be(p, d) 10 Be single 
nucleon transfer reaction leading to the + ground state (lsi/2 neutron transfer) of 
10 Be. We perform zero-range calculations using a modified version of the program 
TWOFNR f2"3fl . The program has been further modified so that the calculated 
adiabatic and quasi-adiabatic radial wave functions, and appropriate zero-range 
strengths calculated for the phase equivalent deep and shallow potentials by Eq. 
(41), can be read in, together with the transferred neutron bound state n ob- 
tained using the superdeep sec/i-squared potential and its phase-equivalent part- 
ner. In the three-body model calculations of the deuteron channel wave function 
(d+ 10 Be), we make use of the global optical potential parameter set of Bechetti 
and Greenlees [24] evaluated at half the final state center of mass energy. The 
spin-orbit interactions are included. The entrance channel proton optical poten- 
tial parameters are taken from [pi]] . The spin-orbit force in the proton channel 
is fixed at 6 MeV. The radial integrals are carried out from to 35 fm in steps 
of 0.1 fm. The maximum number of partial waves used was 30 for both the en- 
trance and exit channels. The spectroscopic factors are set to unity throughout 
the calculations. All calculations presented here are done without non-locality 
corrections. Such corrections for halo transfer are expected to be small because 
they correct the transition amplitude in the nuclear interior, but the long tail of 
the halo wave function makes internal contributions less important. 



5.3 Results and discussion 

It is useful at this stage to remind ourselves again that the deep potential and 
its phase equivalent shallow partner, which are used for calculating the bound 
deuteron and transferred neutron wave functions, are constructed to have iden- 
tical phase shifts so that any difference in the transition amplitudes, and in the 
cross section angular distributions, is attributed entirely to the corresponding 
wavefunctions of the partner potentials. 

The result for the adiabatic model cross section angular distributions for the 
reaction involving the original deep potential description and its comparison with 
that obtained by means of the PEP description is given in Fig. 5. The figure indi- 
cates almost complete coincidence of both curves. This may be understood from 
the following analysis. The two-body supersymmetric partner potential depen- 
dence in the differential cross section calculations, originates from two terms: 
the zero-range constant Dq, and the transition amplitude \M(p,d)\ 2 through 
the deuteron and bound neutron ground state wave functions. The calculated 
transfer strengths, for the deep potential description of the n — p interaction is 
D 2 (deep) = 15792 MeV 2 fm 3 while D 2 (pep) = 15980 MeV 2 fm 3 for the shallow 
partner description. It is also worth stressing that, from the results obtained in 
sections 3 and 4, the transition from the deep potential to the repulsive core in- 
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teraction (PEP) does not significantly alter the outer part of the neutron ground 
state wave functions appear in the transition amplitude. And in addition, the 
neutron halo wave function makes internal contribution less important. Therefore 
M dee P(p, d) « \Mw(p, d)\ 2 , and the ratio of the cross sections 

d<J d p 7 P ... D 2 (deep) _ 15792 MeV 2 fm 3 
díl p J ~ Dl{pep) 15980 MeV 2 fm 3 [ } 

= 0.988 w 1 , 

which leads to the coincidence of the results. 



6 Conclusion 

The properties of the deep nuclear interaction have been investigated by con- 
structing explicitly phase-equivalent potentials freed from the unphysical deeply 
bound states of the former. We have seen that the resulting central potentials 
have to be repulsive and singular at small distance in order to preserve the energy 
behavior of the phase shifts, and they present a shallow attractive part of inter- 
mediate range. Our reconstructed potentials (PEP) have led to relative motion 
wave functions very similar to those generated by the deep potentials outside 
the core region, but which lack the small distance radial node. Both types of 
potentials are therefore expected to display rather different off-shell behaviors, 
and presumably lead to qualitatively different results. However, there is no con- 
siderable discrepancy between the rms calculation results of these quite different 
two-body interaction descriptions has been found. Nevertheless, if experiment 
actually does require a reduction in observables then a short range non-local 
contribution to the potential, like the deep potential, appears to be necessary. 

We have used phase equivalent two-body potentials with a different number of 
bound states considering the 11 Be(p : d) w Be reaction at 35 MeV, and compared 
the calculated corresponding physical observables. Investigation of the conse- 
quences of using these completely phase equivalent two-body potentials for the 
description of weakly bound deuteron and 11 Be nuclei in three-body calculations, 
based on the adiabatic approach has led us to almost indistinguishable results. 
Due to the large spatial extension of halos, involving the simplest halo nucleus 
the deuteron and n Be as a one-neutron halo system, the probability is by def- 
inition very small at small distances. Hence, we conclude that the short range 
behavior of the corresponding wavefunctions for the deep and phase-equivalent 
shallow potentials, which coincide at large distances but differ at small distances 
by the additional node appearing inside the core by use of the deep potential, is 
not significant for the analysis of such reactions. 

In sum, the supersymmetric formulations used through the present calcula- 
tions have dealt in general with the Pauli principle for the weakly bound systems. 
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Assuming that the two-body potentials have Pauli forbidden states, one can then 
construct easily and use the phase equivalent partners without these forbidden 
states. At small distances the lowest levels of the original deep potentials corre- 
spond to identical fermions occupying the same states. Removal of these terms 
therefore forces the partides to occupy higher-lying orbits and thereby introduc- 
ing the necessary repulsion preventing violation of the Pauli principle. In con- 
clusion, this method to exclude the Pauli forbidden states in the weakly bound 
systems has firm mathematical and numerical foundations. It is a practical and 
accurate alternative to the other existing methods, such as the work described 
in Ref. ]25| where an analytical s-wave potential with one bound state, which is 
the most important case in the practical applications for halo states, has been 
introduced. We note that the application of this potential, with the appropriate 
choice of the parameters involved, to the weakly bound deuteron and 11 Be nuclei 
has led to the similar results to those obtained by the shallow supersymmetric 
phase-equivalent partner potentials. 

The valuable discussions with D. Baye, W. Van Dijk, and M. W. Kermode 
are here gratefully acknowledged. 
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Figure Captions 



Fig. 1-a. Superdeep potential Vi(r) (sòlid line) for the deuteron system 
and its supersymmetric partners V^r) (non-PEP, dotted line) and V^r) (PEP, 
dashed line) function of radius r. 

Fig. 1-b. Comparison of the n — p central potential of the superdeep sech- 
squared Vi(r) (sòlid line) potential and its PEP potential V 3 (r) (dashed line) with 
the central Reid Soft Core (dotted line) interaction. 

Fig. 2. The first two eigenstates, n = for the ground state (sòlid line) 
and n — 1 for the first excited state (dotted line) of the original hamiltonian 
with the superdeep two-parameter potential for the deuteron. The wave function 
illustrated by dashed line represents the ground state of the SUSY PEP, V 2 (r). 

Fig. 3. The same as Fig. 1-a, but calculations have been carried out for the 
11 Be system. 

Fig. 4. The 11 Be ground-state radial wave functions \E , i™ _1 ' ) (À, /3, r) (sòlid 
curve , ¥^-°\À,(3,r) (dotted curve and ¥^-°\À,f3,r) (dashed curve) all nor- 
malized to one. 

Fig. 5. Calculated differential cross section angular distributions within 
the adiabatic model for the n Be(p, d) 10 Be (g.s.) reaction at 35 MeV using the 
deep and shallow (PEP) two-body potential descriptions for the weakly bound 
deuteron and u Be nuclei. 
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